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Abstract 

We use a sample of 19 Gamma Ray Bursts (GRBs) that exhibit single-peaked 
optical light curves to test the standard fireball model by investigating the relation- 
ship between the time of the onset of the afterglow and the temporal rising index. 
Our sample includes GRBs and X-ray flashes for which we derive a wide range of 
initial Lorentz factors (40 < T < 450). Using plausible model parameters the typical 
frequency of the forward shock is expected to lie close to the optical band; within 
this low typical frequency framework, we use the optical data to constrain e e and 
show that values derived from the early time light curve properties are consistent 
with published typical values derived from other afterglow studies. We produce 
expected radio light curves by predicting the temporal evolution of the expected 
radio emission from forward and reverse shock components, including synchrotron 
self-absorption effects at early time. Although a number of the GRBs in this sample 
do not have published radio measurements, we demonstrate the effectiveness of this 
method in the case of Swift GRB 090313, for which millimetric and centrimetric 
observations were available, and conclude that future detections of reverse-shock 
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radio flares with new radio facilities such as the EVLA and ALMA will test the low 
frequency model and provide constraints on magnetic models. 

1 Introduction 

With the advent of rapid optical follow-up observations of Gamma Ray Bursts (GRBs) 
(e.g. Mundell et al. 2010, Rykoff et al. 2009), the confirmed lack of bright optical flashes 
from most GRBs challenges a key prediction of the standard fireball model in which a 
reverse shock should produce bright, short-lived optical emission at early time (Meszaros 
& Rees 1999; Sari & Piran 1999; Kobayashi 2000). Although the lack of optical flash 
could be partially due to late observations which are not prompt enough to catch early 
flashes, it is not trivial how to explain events like GRB 090313 which exhibits the onset 
of afterglow without signatures of optical flash. 

At early time, reverse shock emission should dominate optical band and a bright 
optical peak is expected to be observed when a fireball starts to be decelerated. How- 
ever, a distinctive reverse shock component is detected only in a small fraction of GRBs 
(Melandri et al. 2008). Several afterglows show a fattening in the light curves, inter- 
preted as the signature of the rapid fading of reverse shock combined with the gradual 
dominance of forward shock emission (Akerlof et al. 1999; Sari &; Piran 1999). Afterglow 
modeling of such flattening cases implies that the magnetic energy density in a fireball, 
expressed as a fraction of the equipartition value of shock energy, is much larger than 
in the forward shock (but it still suggests a baryonic jet rather than a Poynting-flux 
dominated jet: Fan et al. 2002; Zhang et al. 2003; Kumar & Panaitescu 2003; Gomboc 
et al. 2008). Polarization measurements in a rapid decay phase of GRB 090102 afterglow 
shows the existence of large-scale magnetic fields in the revere shock region (Steele et 
al. 20093). The lack of optical flashes in most GRBs may be due to extreme magnetic 
field properties, either high magnetic energy densities that suppress the reverse shock 
(Gomboc et al. 2008; Mimica et al. 2009) or very low magnetic energy densities that 
cause shock energy to be radiated at higher frequencies than the optical band due to 
synchrotron self-Compton processes (Beloborodov 2005; Kobayashi et al 2007; Zou et 
al 2009). Alternatively the light curve flattening could be the result of refreshed shocks 
and episodes of energy injection (Rees &: Meszaros 1998, Melandri et al. 2009). 

A more conventional model would imply that the reverse shock emits photons at 
frequencies much lower than the optical band. Synchrotron emission is known to be 
sensitive to the properties of emitter. Within this framework, which we term the low- 
frequency model, a single peak in the early time optical light curve is produced when 
both of the typical synchrotron frequencies of forward and reverse shock lie below the 
optical band (Mundell et al. 2007a); the single peak actually consists of photons equally 
contributed from forward and reverse shock, the peak time represents the the deceleration 
of a fireball and hence it provides a direct estimate of the initial Lorentz factor. 

1 Mundell et al. 2007b found no ordered magnetic fields or a very high magnetic energy density in 
the ejecta of GRB 060418. More observations are needed to give a strong conclusion on the nature of 
the ejecta (baryonic versus Poynting flux dominated) and the distribution of magnetization degree. 
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In this paper, we discuss the lack of optical flashes in the context of the low-frequency 
model. GRB 090313 is a typical case of a burst that displays a rising and falling light 
curve, little temporal structure, no strong spectral evolution and well-monitored multi- 
wavelength behavior from early times. Here, we analyse its multi-wavelength properties, 
place it into the wider context of GRBs with single optically peaked light curves and use 
the characteristics of the full sample to test the low-frequency model and its predations 
for radio light curve evolution. Throughout the paper we use the following conventions: 
the power-law flux density is given as F(y,t) oc t~ a v~^, where a is the temporal decay 
index and j3 is the spectral slope; a positive value of a corresponds then to a decrease 
in flux, while a negative value indicates an increasing in time of the observed flux. We 
assume a standard cosmology with Hq = 70 km s _1 Mpc -1 , $7 m = 0.3, and Qa = 0.7; 
and all uncertainties are quoted at the la confidence level (cl), unless stated otherwise. 

2 Observations 

On 2009 March 13 at 09:06:27 UT (=T0) the Burst Alert Telescope (BAT; Barthelmy 
et al. 2005) onboard Swift triggered on GRB 090313 (Mao et al. 2009a). The BAT light 
curves showed a series of multiple peaks with the emission starting before T0-100 s and 
a Tgo in the 15-350 keV band starting at ~ TO-3.9 s for a total duration of 78 ± 19 s 
(Mao et al. 2009b). 

Spectroscopic observations performed with the Gemini South telescope provided a 
redshift of z=3.375 for GRB 090313 (Chornock et al. 2009b), later confirmed by from 
VLT with FORS (Thoene et al. 2009) and X-shooter (de Ugarte Postigo et al. 2010; 
who derive a refined redshift value of 3.3736 +/- 0.0004) observations. The estimated 
redshift for this afterglow confirmed again that the near object reported by Berger (2009) 
is indeed too bright to be the host galaxy of GRB 090313. Most likely this extended 
object is one of the two absorbing systems spectroscopically detected (at redshift z=1.96 
or z=1.80) along the line of sight of GRB 090313 (de Ugarte Postigo et al. 2010). Radio 
observations performed with the AMI Large Array (Pooley 2009abc), the VLA (Frail &: 
Chandra 2009) and the WSRT (van der Horst & Kamble 2009ab) confirmed the detection 
and fading nature of the afterglow. 

This event displayed an average 7-ray fluence of ~ 1.4 x 10~ 6 erg cm~ 2 (Mao et al. 
2009b). The redshift of the burst (correspondent to a luminosity distance of ~ 2.9 x 10 4 
Mpc) resulted in an isotropic energy estimate of ~ 3.4 x 10 52 ergs in the 15-150 keV 
observed bandpass. 

2.1 Swift/XRT and Swift/UVOT data 

Due to Moon distance observing constraints there were no prompt XRT (Burrows et al. 
2005) and UVOT (Roming et al. 2005) observations. Follow-up observations of the BAT 
error circle were possible only after ~ 27 ks, showing a power-law decay in the X-ray 
(Mao &: Margutti 2009) and a possible marginal detection in the UVOT-v and UVOT-b 
filters (Schady et al. 2009, Mao et al. 2009b). 
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2.2 Optical and Infrared data 

The optical afterglow was discovered by the KAIT telescope (Chornock et al. 2009a) 
and later confirmed by the GROND telescope at equatorial coordinates (J2000) R.A. = 
13H3 m 36.21 s ; Dec = +08°05'49.2" (Updike et al. 2009). The 2-m Faulkes Telescope 
North (FTN) observed the optical afterglow of GRB 090313 starting from 168 s after the 
burst (corresponding to 38 s in the rest frame). Observations continued up to several 
weeks after the burst with FTN, the 2-m Liverpool Telescope (LT) and the 2-m Faulkes 
Telescope South (FTS) (see Tabled]). Late time observations were also performed in 
order to better correct the entire data set from the contribution of the nearby object, 
close to the position of the afterglow. This object was found to have a constant flux 
equal to ~ 1% of the peak flux of the optical afterglow, not affecting the shape of the 
light curve at early time. 

The optical afterglow was observed also with the 1.5m telescope at the Observatorio 
de Sierra Nevada (OSN), the 0.8m IAC telescope, the 1.23m telescope at the Calar Alto 
Astronomical Observatory (CAHA) and the 0.5m Mitsume telescope in the optical bands 
(R and I), plus the 2.5m Nordic Optical Telescope (NOT) and the 3.5m CAHA telescope 
in the near infrared bands (J and K). It was then possible to build the light curve for 
all the filters as shown in Figfl] A log of the observations is given in Table [H where 
we report the mid time, integration time, magnitude and fluxes for all our detections at 
different wavelengths. Afterglow detections reported in GCNs are also shown in Figfl] 

The optical data were calibrated using a common set of selected catalogued stars 
present in the field of view. SDSS catalogued stars were used for r' and i' filters, while 
USNO-B1 R2 and / magnitudes have been used for the R and / filters respectively. 
J and K observations were calibrated with respect to the 2MASS catalog. Next, the 
calibrated magnitudes were corrected for the Galactic absorption along the line of sight 
(Eb-v = 0.028 mag; Schlegel et al. 1998); the estimated extinctions in the different 
filters are A R ~ A r > = 0.074 mag, Ai ~ A v = 0.054 mag, Aj = 0.025 mag, A H = 
0.016 mag and Ak = 0.010 mag. Corrected magnitudes were then converted into flux 
densities, F u (mJy), following Fukugita et al. (1996). Results are summarized in Table 

m 

2.3 Radio, mm and sub-mm data 

Continuum observations at 870 (im were carried out using LABOCA bolometer array, 
installed on the Atacama Pathfinder Experiment (APEXEI) telescope. Data were ac- 
quired on 2009 March 17 and 24 during the ESO program 082.F-9850A, under good 
weather conditions (zenith opacity values ranged from 0.24 to 0.33 at 870/im). Obser- 
vations were performed using a spiral raster mapping, providing a fully sampled and 
homogeneously covered map in an area of diameter ~12', centered at the coordinates 
of the optical afterglow of GRB 090313. The total on source integration time of the 

2 This work is partially based on observations with the APEX telescope. APEX is a collaboration be- 
tween the Max-Plank-Institut fiir Radioastronomie, the European Southern Observatory and the Onsala 
Space Observatory. 
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two combined epochs was ~ 4.6 hours. Calibration was performed using observations of 
Saturn as well as CW-Leo, B13134, G10.62, and G5.89 as secondary calibrators. The 
absolute flux calibration uncertainty is estimated to be ~ 11%. The telescope pointing 
was checked every hour, finding an rms pointing accuracy of 1.8". Data were reduced 
using the BoA and MiniCRUSH softwares. Finally, the individual maps were co-added 
and smoothed to a final angular resolution of 27.6". We obtained a 3a detection upper 
limit of 14 mJy for each of the two epochs. 

The radio afterglow of GRB 090313 was successfully detected by the AMI Large 
Array ~ 2.8 days after the burst (Pooley 2009a) and then monitored up to ~ 47 days 
(Pooley 2009bc) as reported in Table [21 After an initial upper limit at ~ 1.7 days (van 
der Horst & Kamble 2009a) a detection was reported also by the Westerbork Synthesis 
Radio Telescope at ~ 7.6 days (WSRT, van der Horst & Kamble 2009b) and by the Very 
Large Array at ~ 5.9 days (VLT, Frail & Chandra 2009). In the mm band the afterglow 
was detected with CARMA about one day (Bock et al. 2009) and then monitored with 
the Plateau de Bure Interferometer (PdBI) up to ~ 20 days after the burst event. The 
radio observations are reported in Table [2] where the original frequency range of the 
observation has been specified. 

3 Results 

3.1 BAT spectral and temporal analysis 

We re-binned the BAT light curve of GRB 090313 with dt bins of 16.384 s in order to 
better appreciate the long faint tail visible up to 500 s after the burst onset. As reported 
also by Mao et al. (2009b), the mask- weighted light curve (shown in Fig. [2]) displays a 
series of multiple peaks extending long after t=Tgo at a much fainter level. The time- 
averaged spectrum is best fitted by a simple power-law model with a photon index of 
1.91 ± 0.29 (Mao et al. 2009b). 

3.2 Optical/X-ray light curve 

Observations performed with the Faulkes North Telescope, beginning ~ 170 s after the 
burst, showed the optical afterglow rising to a maximum at ~ 1 ks (Guidorzi et al. 
2009). The peak was followed by a decay with windings and flares (possibly due to the 
interaction with the circum-burst material or late time central engine activities) . Around 
3 x 10 5 s, the magnitude became constant in each filter, revealing the presence of an 
underlying object at the position of the optical afterglow. This faint (r 1 = 21.6 ± 0.2 
and i' = 21.1 ± 0.2) and apparently extended object is only 2.3" away from the optical 
afterglow as reported by Berger (2009). It was not possible to separate the contributions 
from the two objects in the late-time co-added observations. 

We model the optical light curve with a broken power-law (to fit the peak up to ~ 10 4 
s) plus an additional component to model the bumps visible after ~ 1.4 x 10 4 s and a 
constant flux to model the behavior at late times. The fit to the component representing 
the optical peak at early time gives: a r i se = —1.72 ± 0.41, adecay = 1-25 ± 0.08 and 
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tpcak = 1060.9 ± 153.6 s. For completeness the parameters of the component modeling 
the sharp bump around ~ 10 4 s are: a r ,bump = —83.8 ± 8.4, ad,bump = 3.0 ± 0.8 and 
t P eak,bump = (14.0±0.3) x 10 3 s, (t/dt) peak ~ 1 (x 2 /dof = 769.4/77 ~ 9.9). The high Xr 2 ed 
for the optical fit is clearly driven by the uncertainty of the bump fit and the variability 
of the data around ~ 10 5 s. However this does not affect the goodness of the fit for 
the smooth early time behavior, where the peak (rise and fall) is well constrained with 
negligible variability as shown in Fig. [3j 

Our independent analysis shows that the X-ray light curve of GRB 090313 is well 
fitted by a simple broken power-law with a\ = 0.83 ± 0.49, a<i = 2.56 ± 0.46 and 
^break ~ 9 x 10 4 s (x 2 /dof = 43.17/43 ~ 1.0). The estimated values for a± and oli 
could be the result of flares activity, and the subsequent cessation, in the early XRT 
data. The X-ray light curve and its fit are shown in Fig. [3] together with the composite 
optical/infrared light curve. As we will explain in Section 3.4, the latter has been built 
by re-scaling all the filters with respect to the SDSS i' band. On the bottom panel of 
this figure we show the no-evolution of the optical spectral index flo as derived from the 
fit of the spectral energy distribution. 

3.3 X-ray spectral analysis 

The X-ray spectrum (Fig. HI from the .Sro/i-XRT repository, Evans et al. 2007) can 
be fitted by an absorbed simple power law with a photon index Fx = 2.14jl 'j 4 and an 
absorbing column density Ah = (2.99jZQ]y^) x 10 22 cm -2 , in excess of the Galactic value 
of 2.1 x 10 20 cm" 2 . 

3.4 Spectral energy distribution 

From our data and others published in GCNs we estimate the flux for the infrared (JHK) 
and optical (iV) filters at four different epochs (corresponding to T0+100 s, TO+600 
s, T0+2xl0 3 s and T0+1.6xl0 4 s in the rest frame of the burst). At the redshift of 
the burst (z=3.374) the wavelength of the Lyman-alpha break (121.6 nm) is redshifted 
to 532 nm, that corresponds roughly to the central peak wavelength of the V filter. 
However also the tail of the R filter could be affected by the absorption and for that 
reason we decided to perform the fit of the optical spectral energy distribution only up 
to 2xl0 15 Hz. The results of the fit are shown in Fig. [5] and reported in Table The 
afterglow of GRB 090313 did not display any spectral evolution before and after the 
peak in the light curve. Only a slight and insignificant change of the spectral parameter 
fio is recorded around 3ks (observed frame) after the break. For this reason we built a 
composite optical/infrared light curve (fixing the value of /3o = 1-2) using rigid shifts 
for each filter to report all the fluxes relative to the SDSS-i band. 

4 Discussion 

Here we examine the properties of 19 GRBs including GRB 090313 that exhibit a single- 
peaked optical light curve. Those are all the GRBs with published data that show a 
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clear rise and fall of their optical light curves. The observed and derived properties of 
the sample are given in Table HI In this table we report the parameters of the optical 
peak (a r i se , ctdecayj tpeak and F p ), together with the X-ray decay index (ax) in the 
post optical peak phase □, the duration (T90), redshift (2), initial Lorentz factor T and 
isotropic energy (Ei SO ) for each burst. We have assumed that the optical peak time 
represents the fireball deceleration time. Following equation 1 in Molinari et al. 2007, 
the initial Lorentz factor of GRB 090313 is give by 

r » so „-■/• f y /s (i±i) 3,a ( v^r 3/s (1) 

V3.2 x 10 52 erg/ V4.375; V1060 sj KJ 

where n is the ambient density in protons/cm 3 . For all the bursts in Table H] the ISM 
environment is favored in literature (i.e. Klotz et al 2008, Rykoff et al. 2009, Oates et al. 
2009, Melandri et al. 2009, Greiner et al. 2009); only GRB 080330 is better explained 
by a wind- like medium (Guidorzi et al. 2009). For the wind medium p = AR~ 2 , the 
equation 1 is replaced by T ~ 25 (A/5 x 10 11 g cm- 1 )~ 1 / 4 ( J E/3.2 x 10 52 erg) 1 / 4 [(l + 
z)/4.375] 1 /4( W /i060 s)- 1 /4. 

It is well accepted that the X-ray temporal decay of the majority of GRB afterglow 
can be described by a canonical light curve, where the initial X-ray emission (steep decay) 
is consistent with the tail of the gamma-ray emission, followed by a shallow phase that 
leads into a power-law decay phase (Nousek et al. 2006, Zhang et al. 2006, Tagliaferri 
et al. 2005, O'Brien et al. 2006). In our sample also, no peaks are detected in the X-ray 
light curves, all the X-ray light curves monotonically decay from the beginning of the 
X-ray observations except X-ray flares. It is known that about 50% of GRBs show flaring 
activities on top of the canonical light curve. The narrow structure At/t < 1 indicates 
that it originates from a physically distinct emitting region (e.g. late internal shocks). 
X-ray observations started before an optical peak for GRB 990123, GRB 050730, GRB 
050820A, GRB 060418, GRB 060605, GRB 060607A, GRB 060904B, GRB 070419A, 
GRB 074020, GRB 071031, XRF 080330 and GRB 080810, while it started after an 
optical peak for GRB 061007, GRB 080603A, GRB 080129, GRB 080710 and GRB 
090313. We have no X-ray observations for XRF 020903 and XRF 030418. If an optical 
peak is due to the deceleration of a fireball, X-ray emission from external shocks also 
should peak simultaneously. The tail of the prompt emission or a different emission 
component might mask the X-ray peak. 

Four events: GRB 060418, GRB 060605, GRB 060607A and GRB 060904B show 
X-ray flares around an optical peak, we tested whether the observed optical peaks could 
be explained by the flare emission alone by extrapolating the peak flux of the X-ray 
flare to the optical band assuming a spectral index between the two bands of j3 ~ 1. In 
all cases, the contribution of the X-ray flare to the optical light curve was significantly 
lower than that observed, ruling out a flare origin for the optical peaks. 



3 the value of ax is taken from the literature or from the XRT light curve repository (Evans et al. 
2007). 
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4.1 The origin of the optical peak 

Recent results on the naked eye optical flash from GRB 080319B (Racusin et al. 2008; 
Bloom et al. 2009), where the observed optical peak coincided in time with the prompt 
gamma-ray emission, provided motivation to consider that the prompt gamma-ray emis- 
sion is Inverse Compton (IC) of the optical flash. The dominance of IC cooling could 
lead to the lack of prompt optical flashesQ However, the basic problem of such IC 
model is that if the low-energy seed emission is in the optical, while the observed soft 
gamma-ray spectrum is the first IC component, then second IC scattering would create 
a TeV component. The second IC component in the TeV range should carry much more 
energy than the soft gamma-ray components. This could cause an energy crisis problem, 
possibly violating upper limits from EGRET and Fermi (Piran et al. 2009). 

Rykoff et al. (2004) suggested a model in which single-peaked light curves are caused 
by GRB radiation emerging from a wind medium surrounding a massive progenitor. 
This model suggests that the rise of the afterglow observed in the optical band can be 
ascribed to extinction and the emission can be modeled with an attenuated power-law. 
A consequence of this model is that at very early times some afterglows will rise very 
steeply and the extinction observed in the optical band should be much greater that in 
the infrared band. As shown in FigJH] we see a very steep rise only for GRB 061007, 
however for this burst as for the other bursts on that figure, we do not have data to 
model the peak in the infrared band. If we fit the afterglow peak of GRB 090313 with 
an attenuated power-law function (equation 1 in Rykoff et al. 2004) we find values of 
the decay index and the attenuation time scale (a = 1.15 ± 0.03 and /3 t = 1097 ± 117 
s) consistent with the decay index a obtained in section 3.2. With this /3t we derive a 
mass loss rate (~ 10~ 3 M yr" 1 ) which is slightly higher than what is usually suggested 
for GRB progenitors. The Lorentz factor that we assumed for this estimate is obtained 
from the peak time based on the wind model; T based on the ISM model is higher and 
it would results in a higher mass loss rate. This is a similar result to the one found by 
Rykoff et al. for GRB 030418. As the majority of the GRBs in our sample rise slowly 
or with comparable a r i se with respect GRB 090313 this will imply a higher mass loss 
rate for all those bursts. This model will be further tested with future simultaneous 
optical/IR light curves obtained at early time. 

If the observed peak is due to the passage of the typical frequency of the forward shock 

4 The full discussion on IC cooling effects (e.g. Nakar, Ando & Sari 2009) is beyond the scope of 
this paper. We here give a rough estimate on how much e_e would be necessary to suppress an optical 
flash. We assume that the typical frequency of the reverse shock is in the optical band, and that the 
shock emissions are in the fast cooling regime. If the IC cooling is not important, the luminosities would 
peak at the typical synchrotron frequencies, and the luminosities would be comparable at the onset of 
afterglow. The flux ratio is about F in the optical band (Kobayashi & Zhang 2003) . If the IC cooling is 
the dominant cooling mechanism of the electrons in the shock regions, the bulk of the shock energy is 
radiated in high energy radiation (possibly the 1st scattering component for the forward shock and the 
2nd scattering component for the reverse shock). The optical flux ratio could be reduced roughly by a 
factor of (e e /es) 1/6 (Kobayashi et al. 2007). A very small e_e ~ e e /r B is required to explain the lack of 
optical flashes. In the slow cooling regime, the Compton parameter is smaller for a given ratio e e /es, 
the required eb could be even smaller. 
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through the optical band, we would expect much slower rise (ft r i se ~ -0.5) and strong 
color evolution around the peak. These are not consistent with GRB 090313 observations 
(a r ise ~ -1-7 and no color evolution). If the optical peak is due to the deceleration of 
a fireball, the typical frequency of the forward shock z^ m f s should be below the optical 
band at the onset, otherwise, the forward shock emission slowly rises until the typical 
frequency crosses the optical band. Actually when this condition: ^ m ,fs(ipeak) < ^optical 
is satisfied, the forward and reverse shock emission peak at the same time, and produce 
a single peak (Mundell et al. 2007a). We here consider such a low- frequency model in 
detail. 

The onset of the afterglow is expected to occur immediately after the prompt emission 

if the reverse shock is in the thick shell regime, while there should be a gap between the 

prompt gamma-ray emission and the onset if the reverse shock is in the thin shell regime 

(Sari 1997). At the onset of afterglow, the forward and reverse shock emission rise as 

F cx i 3 and t 3p ~ 3 / 2 , respectively in the thin shell case, while they are as shallower as 
t (3- P )/2 and t i/2 for the thick ghell 

case. If the two emission components are comparable 
at the onset, the rising index could be determined by the shallower component. The 
rising index is expected to be t 3 for the thin shell case, and t 1 ^ 2 or shallower for the thick 
shell case. 

As we will discuss, most optical afterglows are classified into the thin shell case. The 
fireball deceleration time is given by t V eak ~ 90 (1 + z)El f 2 3 n" 1 / 3 r 2 8/3 s where we have 
scaled parameters as £"52 = -Ei so /10 52 ergs and T2 = T/100. At the peak time i pea k, the 
cooling frequency and the typical frequencies of the forward and reverse shock emission 
are given (Sari et al. 1998, Kobayashi &: Zhang 2003) by 

v c ~ 2.6 xlO 18 (l + ^^e^^V^r^Hz, (2) 

~ 5.4 x 10 13 (1 + zrH^KV 172 ^ Hz > (3) 

z, m , rs ~ 5.4 x 10 9 (1 + z)- l el_ 2 e^%n l l 2 Tl Hz, (4) 

where e e -2 = e e /W~ 2 and €3-3 = e#/10~ 3 . For plausible parameters, the typical 
frequency of the forward shock is actually below optical band and the both shock emission 
is in the slow cooling regime. 

The low typical frequencies provide an upper limit to the microscopic parameter e e . 
Requiring that the typical frequency of the forward shock is below the optical band at 
the onset of afterglow, we obtain 

^ iM 1 ' 1 (i +^ ,/4 fef (i^r /4 fef < 5 > 

The estimated values for the upper limit of e e for the GRBs in our sample are reported 
in tabled! The spread of the values of eb is large (from ~ 10 -4 to ~ 10 _1 ) and this could 
still be a significant uncertainty in the upper limits estimates, even if e e do not strongly 
depend from that parameter. GRB 020903 does not constrain e e well, the typical upper 
limit is ~ 0.08, consistent with values from later afterglow modeling (e.g. Panaitescu & 
Kumar 2002). 
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A highly magnetized fireball is another possibility to explain the lack of optical 
flashesll However, Granot et al. (2010) recently argued that in the thin shell case 
the magnetization of the GRB outflow at the deceleration time is not high enough to 
suppress the reverse shock. Most events in our sample are classified into the thin shell 
case. Even if the reverse shock is suppressed by high magnetization, the same condition 
^m,fs(ipcak) < ^optical could be required to avoid slowly rising forward shock emission 
after the onset of afterglow. 

4.2 Reverse and Forward Shocks: Relative Contributions 

In Fig. [6] we plot the light curve rise index (a r i se ) against the time of the peak in the GRB 
rest frame (left panel) and the ratio ipeak/Tgo (right panel). In a recent work Panaitescu 
& Vestrand (2008) classified the optical light curves into 'fast-rising with an early peak' 
and 'slow-rising with a late peak'. In our sample an apparent weak anti-correlation can 
be seen between a r i se and t pea k; however the significance is very low (~ 12%) not allowing 
any firm conclusion about the existence of this anti-correlation. A simple fireball model 
predicts that the dynamics of a fireball is classified into two cases: (1) thin shell fireballs 
(tpeak > T90) produce a sharp peak with rising index a r i se ~ —3; (2) thick shell fireballs 
(^peak ~Tgo) have a wider peak with a r i se ~ — 1/2 (where this value is a limit and the 
rise could be much shallower). As shown in Fig. [6] (right panel), most GRBs in the 
sample are classified into the thin shell case, and the rising indexes are consistent with 
the simple model or shallower. The simple reverse shock model assumes a homogeneous 
fireball. However, as internal shock process requires, the initial fireball could be highly 
irregular. The complex structure of shell or energy injection in the post-prompt phase 
could make the rising index shallower. 

In Fig. El GRB 061007 stands out as a notable exception with Rykoff et al. (2009) 
quoting a peculiarly steep rising index (a r i se ~ -9). Mundell et al. (2007a) showed 
that the afterglow is detected from gamma to optical wavelength, beginning during the 
prompt emission as early as 70 s post-trigger. The softening of the gamma-ray spectral 
index after 70 s further confirms the afterglow onset at this time (Mundell et al. 2007a; 
Rykoff et al. 2009). The gamma ray light curve is dominated by a multi-peaked flare 
between T=20 and 70 s, coincident with the steepest rising part of the optical light curve 
and possible double optical peak. If the optical emission during these prompt gamma-ray 
flares comprises a rising afterglow component with a contemporaneous prompt (flaring) 
component superimposed, the underlying afterglow rising index would be much shallower 
than the observed value. In our small sample, the optical afterglow of XRFs tend to rise 
slowly with a late peak. If we ignore XRFs and the peculiar case of GRB 061007 in fig 
[6l the anti-correlation between the rising index and peak time is very weak or it might 
not exist. 

5 The reverse shock emission might be suppressed for high magnetization: a = B 2 /Anpc 2 ~ 0.1 or 
larger where B and p are the rest-frame magnetic field strength and density, respectively (Mimica et 
al. 2009). Assuming a mildly relativistic reverse shock, the critical magnetization could correspond to 
e B -0.1. 
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GRB 990123 has a clear reverse shock component in early optical afterglow. Our 
low-frequency model is not suitable to discuss this event, because it is considered to 
explain the lack of optical flash. On the other hand, the simple reverse shock model still 
predicts that the rising index is ~ 1/2 for the reverse shock dominant thick shell case. 
We plot GRB 990123 in Fig. [6] also to test the simple model. The discrepancy might be 
due to the irregularity of the fireball. 

An interesting comparison can be done with the peaks detected in the high energy 
band by the Fermi/LAT. Ghisellini et al. (2010) studied the emission observed at ener- 
gies > 0.1 GeV of 11 GRBs detected by the Fermi. They argue that the observed high 
energy flux can be interpreted as afterglow emission shortly following the start of the 
prompt emission. Most events show the onset of afterglow during the prompt gamma- 
ray phase. This is quite a contrast to what we have seen in our sample. The reason for 
this difference might be that Fermi events tend to have very high Lorentz factors, which 
allow to emit high energy photons without pair attenuation, and the events are classified 
into the thick shell. The peak time should be early and comparable to the duration of 
the prompt emission. On the other hand, our sample (early optical observations) might 
be biased towards the thin shell case, because early peaks are technically difficult to 
catch. 



4.3 Radio Afterglow Modeling 

In the low frequency model, the characteristics of an optical peak: the peak time £ pea k 
and the peak flux F p can be used to predict the behavior of early radio afterglow. At 
the onset of afterglow (peak time), the typical frequencies and spectral peaks of reverse 
and forward shock are related as v m ^ s ~ r _2 f mi f s and -F m ax,rs ~ rF maXi f s , respectively 
(Kobayashi & Zhang 2003). Note that F p is a peak in the time domain, while F max 
is a peak in the spectral domain. To produce bright forward shock emission, z^ mi f s 
should be close to optical band and we get v m . rs ~ r~ 2 f opt and F maXiTS ~ TF p . Af- 
ter the original fireball deceleration, the typical frequency and spectral peak behave 
as z^ m rs ^ t-^/ 2 and 

-^max,rs ~ t ■ The typical frequency comes to the radio band 
at t ~ x 2/3 Kpt/Vradio) 2/3 r- 4/3 t P eak ~ 2 x 10 3 r- 4 /% eak and the flux at that time is 
F ~ x~ (3p+1)/6 (^o P t/^adio)- 2/3 r 7 / 3 F p ~ 5 x 10- 4 r 7 / 3 F p where X = *Vf>o P t < 1 is a 
correction factor when v m} f s is well below the optical band and F p = x -fmax,fe> in 
principle, x could be determined from radio observations, T is estimated from the peak 
time tpeak as shown in tabled! In eq (5), the upper limit corresponds to the case of 
X = 1. If X is obtained from radio observations, the right-hand side of the inequality 
with a correction factor of x 1//2 gives the value of e e . 

At low frequencies and early times, self-absorption takes an important role and sig- 
nificantly reduces the flux. A simple estimate of the maximal flux is the emission from 
the black body with the reverse shock temperature (Sari & Piran 1999; Kobayashi & 
Sari 2000). The black body flux at the peak time is 

F v .bb ~ 7r(l + z)v 2 e e m p T \jj^J (6) 
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where R± ~ 2Fci pea k is the observed size of the fireball. This limit initially increases as 
~ t^^, and then steepen as ~ i 5 / 4 after 

^m^s crosses the observation frequency v. In Fig. 
[7] the dashed lines indicate the black body flux limit. Once the reverse shock emission 
becomes dimmer than the limit, the flux decays as ~ i - ( 3p+1 )/ 4 . The combination of 
the increasing limit and decaying flux shapes "radio flare" (Kulkarni et al. 1999). The 
forward shock emission (thin solid) evolves as t 1//2 before the passage of v m ^ s through 
the radio band, and then decays as t~ 3 (?> -1 )/ 4 . The forward shock peak F ~ x ~^ p ~ 1 ^ 2 F p 
should happen around t ~ x 2 ' 3 (^optMadio) 2 ' 3 *peak ~2x 10 3 t pea k s. 

In the case of GRB 090313, assuming x ~ 1> the forward shock peaks in the opti- 
cal band f mi f s ~ 4.6 x 10 14 Hz with a flux density F max f s ~ 2 mJy and a peak time 
corresponding to T ~100; therefore the reverse-shock peak flux at this time occurs at 
^m,rs (tpeak) ~ 46 x 10 9 Hz and is F maXjI . s ~ 200 mJy. Correcting for synchrotron self- 
absorption, results in an observable flux density of ~ 4 mJy after 2.4 hours. After the 
deceleration time the reverse-shock emission in the radio band decays as ~ t~ 2 (dot- 
dashed line Fig. [7J and the emission at 1 day is about ~ 20 /iJy for GRB 090313. 

For GRB 090313, the forward shock emission is expected to peak in the radio band 
around 12 days after the burst (assuming v ra dio = 1-5 x 10 10 Hz), with peak flux of 2 mJy 
(solid line Fig. [7J. Taking all these factors into account, the resultant expected light 
curve of the radio afterglow of GRB 090313 is shown in Fig. The expected 15 GHz 
and 100 GHz light curves (thick lines) are reasonably consistent with the observations. 
The deviation of the 15 GHz estimates from the observations might be partially due 
to a simplified synchrotron spectrum which is described by a broken power law. The 
deviation of the 100 GHz point around 1 day is much more apparent and might be 
due to an additional emission component (e.g. late time central engine activity). Since 
a realistic synchrotron spectrum is rounded at the break frequencies, a more accurate 
estimate should give a light curve rounded at the peak time. However, if this is the case, 
our simple model further underestimates the 100 GHz flux. It is interesting that the 15 
GHz flux decays very slowly up to few ten days while X-ray afterglow displayed a steep 
decay around ~ 1 day, as shown in Section 3.2. This might indicate different origins 
(e.g. emission regions) for the two; the 5a = Q2 — a ra dio > 2 is indeed too large to be 
explained assuming that the cooling frequency lies at the X-ray frequencies at that time. 

In Fig [51 we show radio light curves expected for our sample, which are evaluated 
by using early optical observations. GRB 990123, XRF 020903, GRB 030418, GRB 
060607A, GRB 070420 and GRB 080810 are excluded in the radio afterglow estimates. 
Since GRB 990123 clearly shows a reverse shock component in the early afterglow, it is 
not consistent with our model assumption. For the other five events, the optical peak 
time or peak flux was not well constrained. In future, we should be able to estimate 
radio afterglow light curves in real time as soon as a single peaked optical light curve is 
detected. Depending on the Lorentz factor at the time of the peak and on the energetics 
of the burst, the shape of the radio will slightly change, displaying an early peak/flash 
at ~ 0.1 days and later on the peak of the forward shock in the radio band peaking 
at about 2-10 days after the burst. Diffractive scintillation might make the detection 
of radio flares difficult if the amplitude of flares are order of unit. In cases similar 
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to GRB 061007, in which the optical forward and radio reverse shocks peak at early 
time and the forward shock flux is large, the radio peak due to the the passage of the 
forward shock typical frequency is expected to be very bright. Liang et al. (2010) 

2/7 i 

suggest a correlation, such that T oc E-J Q , therefore Fig. [8] can also be viewed in terms 
of increasing E; so . The scatter in this correlation, however results in an over-prediction 
of the the initial Lorentz factor for GRB 090313 (T ~ 130) compared with the value 
calculated directly from the light curves. 

5 Conclusions 

We have analysed multiwavelength observations of GRB 090313 and similar 18 GRBs 
which exhibit a single-peaked optical light curve. We have compared prompt and af- 
terglow properties to test the standard fireball model with amended microphysics pa- 
rameters. The goal of the study was to understand the origin of single optical peaks 
in afterglow light curves and to explain the surprising lack of bright optical flashes 
from reverse shocks that were predicted from the standard fireball model. Within this 
amended standard model, which we term the low-frequency model, a single peak in the 
early time optical light curve is produced when the typical synchrotron frequencies of 
shock emission lie below the optical band. We have shown that this condition is satis- 
fied with plausible microphysics parameter e e ; the single peak consists of forward and 
reverse shock emission components, the peak time represents the initial deceleration of 
the fireball at the onset of the afterglow and the reverse shock emits most photons at 
frequencies below the optical band. We find that: 

• In the case of GRB 090313, no spectral evolution was observed at the time of the 
optical peak, the peak is considered to represent the onset of the GRB afterglow 
(or fireball deceleration) and the initial Lorentz factor of the ejecta was derived 
r ~ 80. The Lorentz factors that were similarly derived for the other GRBs and 
XRFs in the sample cover a wide range 40 < T < 450. 

• The rising indexes of most optical light-curves are consistent or shallower than 
the value of F ~ i 3 expected in the standard model. Although a simple reverse 
shock model assumes a homogeneous fireball, the internal shock model requires a 
highly irregular fireball. At the end of the prompt gamma-ray phase, the fireball 
might still have an irregular structure. The irregularity in the density distribution 
or energy injection in post-prompt gamma-ray phase could make the rising index 
shallower than the expected value. In the small sample, the optical afterglow of 
XRFs tend to rise slowly with a late peak. 

• We constrained the value of e e for the single-peak events, found an average value 
of < 0.08 for the whole sample. The values derived from early time light curve 
properties are consistent with published values derived from late-time afterglow 
modeling. However the large spread of values for cb could affect the estimates of 
the upper limit for e e . 
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• Using the observed optical properties for our sample of GRBs, we predicted the ra- 
dio afterglow light curves for the low-frequency model. Synchrotron self-absorption 
is important at early times in shaping the radio light curve and masking the reverse 
shock emission. This could result in an early detectable peak around ~ 0.1 days, 
though prompt radio observations might be challenging. The forward shock peaks 
later around 2-10 days after the burst. It is important to note that high energies 
and Lorentz factors (as in the case of GRB 061007) could produce bright optical 
and radio afterglows. We demonstrate the effectiveness of this method in the case 
of GRB 090313. This is important for future observations of GRBs afterglow in 
the radio band with new facilities such as the EVLA, ALMA and LOFAR. The 
latter will have a very large field of view, and prompt radio observations could be 
possible. However, LOFAR will operate at low frequencies (below 250 MHz) and 
since synchrotron self-absorption limit F^bb ~ 1,2 is much lower, it could be still 
difficult to catch prompt optical flares. Current radio sensitivities of 50 /xJy are 
already adequate for detecting reverse and forward shock peaks but with predicted 
sensitivities as low as 2.3 /iJy in a 2-hour integration (Chandra et al. 2010) all 
radio light curves in our sample would be easily observed from early to late time 
with instruments such as the EVLA and ALMA. 
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Figure 1: Multi-band light curve in the observed frame for GRB 090313. Filled symbols 
are out data while GCN data (open symbols) are from: Klotz et al. 2009 (TAROT), 
Nissinen 2009 (Gras-04), Vaalsta & Coward 2009 (Zadko), Perley 2009 and Perley at 
al. 2009 (Gemini), Cobb 2009 (SMARTS), Updike et al. 2009 (GROND), Morgan et al. 
2009 (PAIRITEL). 



19 



0.06 
0.05 
0.04 
0.03 
0.02 
0.01 

o —ME 

-0.01 
-0.02 

-100 100 200 300 400 500 
Time from BAT trigger [s] 

Figure 2: Gamma-rays 16 s binning light curve for GRB 090313 as observed by BAT. 



D 

'tn 
u 



1 



20 



3 



100 
10 

1 

0.1 
0.01 
0.001 

0.0001 

-p 

•H 1.0 

o 1.6 

<o 

*j 1-3 
pa 

1.0 



1 5 GHz 

'comp 
f X*100 

Fit 



; x i 



f, 



const 



4- + i » 



+ 



+ 



100 1000 10000 100000 1e+06 1e+07 1e+08 1e+09 

t since GRB trigger [s] 



Figure 3: Composite optical/infrared light curve (blue symbols) for the afterglow of GRB 
090313. The X-ray light curve (black symbols) is well fitted by a broken power-law. See 
the text for the details about the fit of the composite optical/infrared light curve. In 
the lower panel the evolution of the optical spectral index is shown. 
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Swift-XRT PC spectrum of GRB 090313 



■ i 1 1 1 1 — rp - 1 1 1 


— i 1 — i — i — ■ 

— i 1 — i ~p"j 


= — i — i — i — i — i — \— i 1 hj — 


^ 1 ■■ 


. J- t . . + 





• 1 I 1 I I ±_l 1 ■ T . — 1-. 

0.5 1 2 5 

Energy (keV) 



Figure 4: The X-ray spectrum in the time interval [TO+26779, TO+46710] s is well fitted 
by a simple absorbed power law (Evans et al. 2007). 
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Figure 5: Multi epoch spectral energy distribution. Times are in the rest frame of the 
GRB. 
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Figure 6: Left: rise index (a r ise) vs time of the peak in the GRB rest frame. Right: rise 
index (a r i se ) vs tpeak/Tgo- Reverse shocks can be classified into two classes: 1) t pea k > 
T90, thin shell case, sharp rise; 2) t pea k ~ T90, thick shell case, slow rise. The vertical 
line shows where t pe ak = T90. Dashed horizontal lines represent the asymptotic values 
for the two cases (see text for details). The absolute values of a r ise are plotted. 
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Figure 7: Expected light curves at 100 GHz (left) and 15 GHz (right). The reverse shock 
evolution is shown as a dot-dashed line, the self-absorption curve as a dashed line and 
the forward shock evolution as a solid thin line. The thick line represent the expected 
light curve for GRB 090313. Circle and square points represent observed data in the 
mm (100 GHz) and radio (15 GHz) band respectively. 
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Figure 8: Schematic showing general dependencies for different components (a); pre- 
dicted light curves in the radio band, u = ^vla = 8.5 GHz (b) for the sample GRBs, 
derived from their observed optical properties listed in Tabled! The brightest blue curve 
represent the peculiar high energetic case of GRB 061007 for which a very bright radio 
afterglow is expected, see text for more details. 
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Table 1: Observed magnitudes and fluxes for the optical/infrared afterglow of GRB 090313. At is the delay since the burst 
event in the observer frame. Magnitudes are not corrected for Galactic absorption while F v are the absorption corrected 
converted flux densities. 
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FTN 


14.71 


60.0 


Rc 


16.07 ±0.10 


1.751 ±0.134 


FTN 


16.71 


60.0 


i> 


15.31 ± 0.08 


2.922 ±0.129 


FTN 


21.88 


120.0 


Rc 


15.67 ±0.10 


1.897 ±0.193 


FTN 


24.95 


120.0 


i' 


15.40 ±0.05 


2.689 ±0.129 


FTN 


32.21 


180.0 


Rc 


16.02 ±0.11 


1.419 ±0.154 


FTN 


26.20 


180.0 


i' 


15.68 ± 0.05 


2.078 ±0.129 


FTN 


89.28 


30.0 


Rc 


17.62 ±0.20 


0.348 ± 0.065 


FTN 


93.04 


10.0 


i> 


17.03 ± 0.17 


0.599 ±0.129 


FTN 


95.70 


30.0 


Rc 


17.47 ±0.17 


0.399 ± 0.063 


FTN 


97.10 


30.0 


i> 


16.79 ± 0.10 


0.747 ±0.129 


FTN 


100.54 


60.0 


Rc 


16.92 ±0.13 


0.663 ± 0.079 


FTN 


102.44 


60.0 


i' 


16.65 ± 0.08 


0.850 ±0.129 


FTN 


107.70 


120.0 


Rc 


17.51 ±0.11 


0.385 ± 0.039 


FTN 


110.56 


120.0 


i' 


17.09 ± 0.09 


0.567 ±0.129 


FTN 


117.78 


180.0 


Rc 


17.78 ±0.14 


0.300 ± 0.039 


FTN 


121.75 


180.0 


i' 


17.43 ± 0.13 


0.415 ±0.129 


FTN 


127.19 


120.0 


Rc 


18.20 ±0.22 


0.204 ±0.042 


FTN 


130.82 


120.0 


i' 


17.10 ±0.11 


0.562 ±0.129 


FTN 


137.99 


180.0 


Rc 


17.99 ±0.18 


0.247 ± 0.041 


FTN 


141.99 


180.0 


i' 


17.18 ± 0.10 


0.522 ±0.129 


FTN 


204.80 


300.0 


Rc 


18.14 ±0.10 


0.216 ±0.020 


FTN 


169.90 


900.0 


i' 


17.21 ± 0.10 


0.508 ±0.129 


FTN 


210.18 


300.0 


Rc 


18.23 ±0.10 


0.198 ±0.018 


FTN 


186.03 


900.0 


i> 


17.78 ± 0.13 


0.300 ±0.129 


FTN 


215.56 


300.0 


Rc 


18.49 ±0.12 


0.156 ±0.017 


FTN 


336.59 


300.0 


i' 


17.21 ± 0.15 


0.508 ±0.129 


FTN 


220.93 


300.0 


Rc 


18.43 ±0.10 


0.165 ±0.015 


FTN 


352.72 


1500.0 


i' 


17.26 ± 0.16 


0.485 ±0.129 


FTN 


226.31 


300.0 


Rc 


18.37 ±0.12 


0.174 ±0.019 


FTN 


86288.09 


600.0 


i' 


20.94 ±0.22 


0.0163 ± 0.0034 


FTN 


231.69 


300.0 


Rc 


17.53 ±0.15 


0.378 ± 0.052 


FTN 


132306.85 


1800.0 


i' 


21.12 ± 0.10 


0.0136 ± 0.0012 


FTN 


306.69 


900.0 


Rc 


17.40 ±0.10 


0.426 ± 0.039 


FTN 


865.36 


900.0 


i' 


19.01 ± 0.11 


0.097 ±0.010 


LT 


322.82 


900.0 


Rc 


17.27 ±0.15 


0.480 ± 0.066 


FTN 


880.96 


900.0 


i' 


18.73 ± 0.08 


0.125 ±0.009 


LT 


1407.83 


1800.0 


Rc 


19.85 ±0.11 


0.0446 ± 0.005 


FTN 


1010.82 


1800.0 


i' 


19.29 ± 0.25 


0.074 ±0.017 


LT 


86255.28 


600.0 


Rc 


21.58 ±0.15 


0.0091 ±0.0013 


FTN 


1156.30 


900.0 


i' 


18.86 ± 0.05 


0.111 ±0.005 


LT 


86265.98 


600.0 


Rc 


21.63 ±0.13 


0.0086 ± 0.0010 


FTN 


1169.30 


900.0 


i' 


18.84 ±0.07 


0.113 ±0.007 


LT 


86277.15 


600.0 


Rc 


21.75 ±0.15 


0.0077 ± 0.0010 


FTN 


2431.86 


7200.0 


i' 


20.25 ± 0.07 


0.0309 ± 0.0020 


LT 


131244.99 


1200.0 


Rc 


21.40 ±0.64 


0.0107 ± 0.0067 


FTN 


3870.60 


7200.0 


i' 


20.43 ± 0.08 


0.0261 ± 0.0020 


LT 


828.90 


1800.0 


r' 


19.73 ±0.10 


0.0498 ± 0.0046 


LT 


27001.76 


600.0 


i> 


21.09 ± 0.15 


0.0142 ± 0.0020 


LT 


983.50 


3600.0 


r' 


18.65 ±0.11 


0.1348 ±0.0137 


LT 


4872.80 


1800.0 


i> 


21.23 ± 0.18 


0.0125 ± 0.0020 


FTS 


1200.25 


1800.0 


r' 


19.60 ±0.23 


0.0562 ± 0.0120 


LT 


6285.20 


1800.0 


i' 


21.10 ±0.20 


0.0138 ±0.025 


FTS 


2507.32 


7200.0 


r' 


20.76 ±0.11 


0.0193 ±0.0019 


LT 


811.10 


900.0 


I 


18.75 ± 0.17 


0.123 ±0.020 


OSN 


3933.26 


7200.0 


r' 


21.47 ±0.10 


0.0100 ±0.0010 


LT 


829.13 


900.0 


I 


18.62 ± 0.10 


0.139 ±0.013 


OSN 


26990.86 


600.0 


r' 


21.62 ±0.21 


0.0087 ±0.0017 


LT 


844.51 


900.0 


I 


18.64 ±0.10 


0.136 ±0.012 


OSN 


4839.42 


1800.0 


Rc 


21.64 ±0.20 


0.0086 ± 0.0016 


FTS 


855.73 


900.0 


I 


18.88 ± 0.12 


0.109 ±0.012 


OSN 


6251.83 


1800.0 


Rc 


21.44 ±0.15 


0.0103 ±0.0014 


FTS 


868.40 


800.0 


I 


18.32 ± 0.15 


0.182 ±0.025 


OSN 


935.04 


600.0 


Rc 


19.31 ±0.12 


0.0734 ± 0.0081 


IAC80 


1203.40 


400.0 


I 


19.23 ± 0.15 


0.079 ±0.011 


OSN 


1167.45 


1380.0 


Rc 


19.74 ±0.26 


0.0494 ±0.0119 


IAC80 


2239.18 


1000.0 


I 


20.15 ± 0.15 


0.033 ±0.005 


OSN 


9695.8 


13800.0 


Rj 


> 23.7 




1.23Caha 


1885.18 


10260.0 


I 


19.77 ±0.26 


0.048 ±0.011 


Mitsume 


11171.1 


9600.0 


R } 


> 23.5 




1.23Caha 


1173.80 


540.0 


J 


17.45 ± 0.20 


0.1708 ±0.0316 


NOT 


1152.38 


1080.0 


K 


15.33 ±0.25 


0.495 ±0.115 


NOT 


8300.50 


3600.0 


J 


19.27 ±0.12 


0.0319 ±0.0035 


3.5Caha 



Table 2: Detected fluxes for the afterglow of GRB 090313 in the radio and mm bands. 



At 


Frequency Range 


F v 


Ref. 


(days) 


(GHz) 


(mjy) 




1.01 


92.5 


4.00 ±0.60 


CARMA, GCN 9005 


1.71 


4.9 


0.026 ± 0.038 


WSRT, GCN 9000 


2.75 


14.5-17.5 


0.800 ± 0.080 


AMI 


3.68 


14.5-17.5 


0.882 ±0.077 


AMI 


4.57 


105.38 


1.681 ±0.153 


PdBI 


4.67 


14.5-17.5 


0.815 ±0.129 


AMI 


5.63 


228.00 


0.605 ± 0.507 


PdBI 


5.85 


8.46 


0.269 ±0.031 


VLA, GCN 9011 


6.78 


14.5-17.5 


0.718 ±0.097 


AMI 


7.40 


4.9 


0.165 ±0.030 


WSRT, GCN 9016 


9.67 


14.5-17.5 


0.655 ± 0.069 


AMI 


13.59 


87.205 


0.666 ±0.126 


PdBI 


16.74 


14.5-17.5 


0.569 ±0.140 


AMI 


19.54 


110.00 


-0.206 ±0.304 


PdBI 


36.66 


14.5-17.5 


0.427 ±0.100 


AMI 


46.63 


14.5-17.5 


0.080 ±0.130 


AMI 
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Table 3: Fit results of the spectral energy distributions of GRB 090313. t r f is the time 
of the SED in the rest frame while t b s is the corresponding time in the observed frame. 



SED 


t rf 


*obs 


Po 




(s) 


(s) 




1 


100 


437.5 


1.19 ±0.02 


2 


600 


2625 


1.21 ±0.05 


3 


2000 


8750 


1.39 ±0.07 


4 


1.6 x 10 4 


7 x 10 4 


1.35 ±0.30 
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Table 4: GRBs with detected optical peaks. All the redshifts are spectroscopically confirmed but the redshift for GRB 
070420, for which the pseudo-z from [15] and the photo-z from [6] are reported. References: [1] Akerlof et al. 1999; [2] 
Galama et al. 1999; [3] Bersier et al. 2006; [4] Sakamoto et al. 2004; [5] Rykoff et al. 2004; [6] Oates et al. 2009; [7] Perri 
et al. 2007; [8] Cenko et al. 2006; [9] Vestrand et al. 2006; [10] Molinari et al. 2007; [11] Rykoff et al. 2009; [12] Mundell et 
al. 2007a; [13] Schady et al. 2007; [14] Melandri et al. 2009; [15] Klotz et al. 2008; [16] Kruhler et al. 2009a; [17] Greiner et 
al. 2009; [18] Guidorzi et al. 2009; [19] Guidorzi et al. 2009b; [20] Falcone et al. 2006; [21] Page et al. 2006; [22] Stratta et 
al. 2009; [23] Page et al. 2009; [24] Golenetskii et al. 2007; [25] Kruhler et al. 2009b. Where possible, the value of E iso has 



00 

o 



been taken from Amati et al. (2008) and Rossi et al. (2008) 


GRB 


"rise 


"decay 


tpcak 

00 




A V 
mag 


Too 

00 


z 


r 


ee 

Upper Limit 


(mjy) 


E i so 
(10 52 erg) 


Ref. 


990123 


> -3.5 


2.1 ± 0.2 


50 ± 15 


1.46 ± 0.04 


0.053 


63.3 ± 0.3 


1.61 


S3 450 


0.004 




229 ± 37 


[lj,[2j 


020903 
030418 


-4.6 ± 3.7 
-0.62 ± 0.13 


1.8 ± 3.0 
1.34 ± 0.06 


~ 10 5 
2401 ± 303 




0.098 
0.077 


9.8 ± 0.6 
135 ± 5 


0.25 

< 5 




18.466 


0.10 ± 0.05 
0.40 ± 0.10 


(2.4 ± 0.6) X 10 -3 


[3], [4] 
[5] 


050730 


0.15 ± 0.50 


0.89 ± 0.05 


~ 750 


44+ - 14 
-0.08 


0.155 


155 ± 20 


3.97 


a 110 


0.043 


0.57 ± 0.10 


~ 8 


[6], [7] 


050820A 


-0.35 ± 0.02 


0.97 ± 0.01 


~ 600 


0.93 ± 0.03 


0.137 


26 ± 2 


2.615 


S3 145 


0.021 


4.05 ± 0.15 


97.4 ± 7.8 


[8], [9] 


060418 




1.28 ± 0.05 


153 ± 10 


1.4 ± 0.1 


0.702 


52 ± 1 


1.489 


S3 165 


0.014 


8.00 ± 0.50 


13 ± 3 


[10], [20] 


060605 


-0.84 ± 0.12 


1.16 ± 0.06 


475 ± 53 


0.34 ± 0.08 


0.159 


79.1 ± 3.0 


3.78 


~ 110 


0.042 


2.00 ± 0.10 


2 5+ 3 1 


[11] 


060607A 


-3 6+° 8 
° -1.1 


1 27+ - 16 
1 " 4 '-0.11 


180+;! 


1.09 ± 0.04 


0.089 


100 ± 5 


3.082 


S3 180 


0.014 
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[10], [21] 


060904B 
061007 


_1 7+20 
L - ' -0-7 
-0.82 ± 0.15 
9.0 


1 H + 0' 14 

1.00 ± 0.18 
1.7 ± 0.1 


479. 4+ 2 * 
~ 550.0 
70 ± 15 


0.76 ± 0.04 
0.86 ± 0.05 
1.68 ± 0.02 


0.549 
0.024 


172 ± 5 
75.3 ± 5.0 


0.703 
1.261 


~ 60 
~ 265 


0.093 
0.005 


0.58 ± 0.08 
280.0 ± 10.0 


30+° 19 
86 ± 9 


[11] 
[15] 
[12], [13], [11] 


070419A 


-1.56 ± 0.70 


0.61 ± 0.09 


450 ± 20 


i 97 +0.18 


0.087 


115.6 ± 5.0 


0.97 


~ 60 


0.100 


0.15 ± 0.05 


0.16 


[14] 


070420 


-1.26 ± 0.19 
-0.73 ± 0.14 


0.88 ± 0.09 
1.67 ± 0.15 


~ 200 ± 10 
~ 350.0 


0.23 ± 0.05 


1.561 


76.5 ± 4.0 


1.56 ± 0.35 
3 01 +0.96 
d - U1 -0.68 


S3 230 








[15], [24] 
[6] 


071031 
080129 


-0.7 ± 0.1 
-1.35 ± 0.15 


0.97 ± 0.06 
0.5 ± 0.1 


~ 1000 
~ 12000 


0.99 ± 0.12 
1.5 ± 0.1 


0.036 
3.046 


180 ± 10 

48 ± 2 


2.692 
4.394 


~ 70 
S3 45 


0.102 
0.355 


0.50 ± 0.05 
0.25 ± 0.10 


~ 1.5 
~ 7 


[16] 
[17], [22] 


080330 


-0.4 ± 0.2 


2.02+°;«J 


~ 600 


0.26 ± 0.10 


0.051 


61 ± 9 


1.51 


< 60 


0.109 


0.35 ± 0.05 


0.21 ± 0.05 


[18] 


080603A 
080710 
080810 
090313 


-3.2 ± 0.3 
-1.11 ± 0.03 
-1.32 ± 0.11 
-1.72 ± 0.41 


1.2 ± 0.1 
0.63 ± 0.02 
1.22 ± 0.09 
1.25 ± 0.08 


~ 1600 
1829 ± 19 

~ 100 
1061 ± 154 


1.0 ± 0.1 
1.0 ± 0.2 
1.81 ± 0.20 
0.83 ± 0.49 


0.138 
0.230 
0.087 
0.090 


180 ± 10 
120 ± 17 
108 ± 5 
78 ± 19 


1.688 
0.845 
3.355 
3.375 


~ 55 
~ 40 
~ 260 
~ 80 


0.137 
0.213 
0.007 
0.074 


0.18 ± 0.05 
0.90 ± 0.10 
0.003 ± 0.001 
1.80 ± 0.20 


1.5 ± 0.5 
~ 0.56 

~ 30 
~ 3.2 


[19] 
[25] 
[23] 
This work 



